function [n_nodes,epsi_nodes,weight_nodes] = GaussHermiteNdim(GH_points,ne,stdShocks)

% The grid for numerical itegration 
[GH_e, GH_w] = GaussHermite(GH_points);
tmp_e        = cell(1,ne); 
tmp_w        = cell(1,ne); 
for i=1:ne 
    tmp_e(i) = mat2cell(GH_e*stdShocks(i),length(GH_e),1); 
    tmp_w(i) = mat2cell(GH_w,length(GH_w),1); 
end 
epsi_nodes = cartprod(tmp_e); 
weight_nodes = prod(cartprod(tmp_w),2); 
n_nodes    = length(weight_nodes);

end